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ABSTRACT 


This paper discusses optimal searches for a fixed object and the 
rigorous analytical results of discrete search theory are presented. 

They show that the totally optimal, the uniformly optimal, the locally 
optimal, and the fastest searches are identical under not too restrictive 
assumptions. The mathematical formalism is illustrated by an Earth- 
approaching asteroid search and optimal searches for such objects are 
explicitly constructed. The approximation that Earth-approaching asteroids 
are fixed is equivalent to having a very high (> 100 square degrees/hour) 
search rate. Generalizations to other types of astronomical search are briefly 
mentioned. 
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1 . CONTEXT 


Searches for asteroids, especially Earth-approaching asteroids, 

if 

are routinely carried out by an M.I.T. group (Taff and Sorvarl 1980, 

Taff 1980a, b, 1981, various Minor Planet Circulars), Shoemaker and 
Helln (1979), and others. Although techniques differ (video signal 
processing for us versus the traditional procedures utilizing photo- 
graphic plates for other groups), limiting magnitudes differ, and 
search rates differ, all groups have constrained their searches. In 
particular we all tend to look near the opposition point, during the 
New Moon phase, and especially In the winter months. These limitations 
Increase the brightness of the sought after minor planet and decrease 
that of the night sky background both by limiting scattered light and by 
minimizing background sources of light. The questions addressed In this 
paper are "Are these searches optimal? Xi there an optimal search? In 
what sense is It optimal? How can It be executed?" The answers are "No. 
Yes. Several (and they lead to the same search plan). Simply". 

The branch of mathematics that deals with search theory Is operations 
research and I shall assume that the reader is not well -versed in such 
matters. Hence, a large part of this paper Is, necessarily, an introduction 
to the theory of search. There exists an excellent reference on the subject 
by L. D. Stone (1975). In order to ease the transition for the reader from 
Stone's book to this paper I have followed his notation. Proofs and 
supplementary material can be found therein. 

*Myself, D. E. Beatty, R. L. Irelan, R. C. Ramsey, and J. M. Sorvarl. 
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Below I have formulated what is known as the search with discrete 
effort for a fixed target. Asteroid searches are used as a model to 
illustrate the mathematical formalism. Relatively little simplification 
of the physics or astronomy is necessary to do this. Next the results of 
optimal search theory are stated and the optimal search problem is solved. 
Following this optimal searches for Earth-approaching asteroids are 
e_xp1icitly constructed and exhibited. Lastly, generalizations of optimal 
searches in this and in other astronomical contexts are briefly considered. 

2. FORMULATION 

One looks for asteroids on the celestial sphere. In the largest sense 
this fonns the two dimensional search space of the problem. In practice 
we delineate a limited area of the celestial sphere (say above altitude 
30®) that we shall actually search in. Denote this search space by J. 

One searches using a telescope with a finite field of view. In practice 
we always examine an entire field, never a fraction of a field nor more 
than one field at a time. Hence the search space J consists of a discrete 
set of fields of view. Number these by the index j = 1, 2, . . . In 
particular, since the celestial sphere encompasses 47 t steradians, max (j) < <». 

Before the minor planet is found one assigns an a priori target distri- 
bution on the search space J, p: J -> [0,1] (the notation means that p is a 
function defined over the set J which maps elements of J into the domain 
zero to unity inclusive). The target distribution is the a priori proba- 
bility of finding an asteroid in field of view j e J before one starts the 
search. For main belt asteroids a reasonable model for p is p is uniform 
over all heliocentric ecliptic longitudes and over the heliocentric ecliptic 
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latitude range < 10° (or 5° or 20°). For Earth-approaching minor planets, 
both because of parallax effects and the inherent spread of Earth-approaching 
asteroids over orbital element space (particularly in inclination), a reason- 
able model for p is that p is uniform over the topocentric celestial sphere. 

In any case 

E 

jeJ p(j)<l 

When one examines a field of view (whether a video frame or an exposed 
plate) for an asteroid one expends a certain amount of effort trying to 
detect the asteroid. In the photographic case one is looking for the streak 
that the moving minor planet has left. In the video mode one looks for two 
displaced dots from frames taken at different times (and after the stars have 
been electronically subtracted). One may look at the same field of view 
several times. The cost of performing k inspections in the j'th field of 
view is measured by a cost function c(.i,k); J x {0,1,2,...} ^ [0,~]. 

Clearly c(j,0) = 0 V j e J (no effort implies no cost). One could measure 

cost by the time spent examining a field of view plus the time spent in 
moving to the next field of view (this makes c non-local and is not desirable). 
Operationally we always spend the same time in each field of view (more or 
less). Also, because [area (J)]^'^^/slew speed « time spent examining a 
field of view, the non-local element of c is both unimportant and varies little. 
We typically spend 45 seconds examining a field, the telescope slews at 4°/sec, 
and we rarely search more than 500 square degrees per night. Hence [area 
(J) ' ] slew speed = 5.6 seconds. (For photographic searches it's a good 
approximation too because large plates we usually used; i.e., 6°x6°.) 

Thus I shall measure cost by time and, in the instance of the asteroid search, 
specialize to the case when the incremental cost of the k'th examination in 
field of view number j, viz. 
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Y(J.k) " c(j.k) - c(j,k-l) 


is a constant independent of both j (i.e., the telescope is fast or the 
plates are large and all fields of view a.-e treated equally) and k (e.g.» 
the same field of view is equally well inspected each time). 

When one does examine a field of view of the search space looking for 
an asteroid then there is a conditional probability of detecting it on or 
before the k'th inspection of that field of view (given that it is there). 

This function, for field of view number j and examination k, is denoted by 
b(j,k): J X {0,1,2,...} ^ [0,1]. Naturally b(j,0) = OV j e J (you can't 
find it if you don't look for it). From the detection function b one can 
construct the probability of failing to detect the asteroid on the first 
k-1 scrutinizations of field of view number j and then succeeding on the k'th 
one (given that the asteroid is in field of view number j); viz. 

e(j.k) = b(j,k) - b(j,k-l) 

There is a lot of physics and mathematics subsumed in the detection 
function. Clearly it depends on the asteroid's apparent magnitude, the 
background st^r density, the night sky background brightness, the resolution 
element size ot the detector(s), the false alarm probability one is willing to 
accept, how tired one is, etc. Since the ecliptic is unchanging, atmospheric 
extinction can be computed, the Moon's position is known, etc. this is a 
computable function. Indeed we are developing software to realistically do 
so in a physically correct way. Operationally, for a fixed set of external 
parameters, our detection probability has the shape shown in Fig. 1 where 
m^ is our quoted limiting magnitude (e.g., where the probability of detection is 
50%). The form shown in the diagram will be used to compute the ''•'timal search 
plans given below. 


7 



Finally I need to define a search plan . A discrete search plan is 
a sequence C ■= (?i» ? 2 » ^ 3 .***) which tells the searcher to first look in 
cell C.j; if the asteroid is found there then terminate the search but if 
the minor planet isn't found then look next in field of view ^ 

global way to describe this is by a function which specifies the allocation 
of effort devoted to each field of view j. To this end define f: J ■+ [0,“); 
f(j) is the r umber of examinations in field of view number j. 

Above I referred to searches for a fixed target. Clearly the minor 
planets we are trying to find are moving. Earth-approaching asteroids 
ca»^ have geocentric angular speeds of a degree per day (or more). I've 
made the assumption that even these objects are fixed when compared to our 
search rate. The mathematical formulation of this app»'cximation is [area 
(J)/search rate]* asteroid angular speed « field of view. For our parameters 
[area (J) = 500 square degrees, search rate = 100 square degrees/hour, 
asteroid angular speed = l®/day, field of view = 2°] an asteroid could 
traverse only one twelfth of a field of view before the entire search is 
completed. Hence the real problem fits into the formalism reasonably well. 
For faster moving minor planets the optimal search plans developed below need 
to be corrected for the asteroid's motion. 

3. OPTIMAL SEARCHES 

Given the cost of searching field of view number j a total of k 
times, c(j,k), the total cost of performing the search plan C with allo- 
cation f is 

E 

CLf] = jeJ c(j,f(j)) 

E 

The total number of examinations over all fields of view is jeJ f(j). 
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similarly the total probability of minor planet detection with this 
allocation of effort is P[f], 

P[f] - .M P(j)b(j.f(j)) 

where b(j,k) is the conditional probability of finding the asteroid in 
field of view number j after k examinations of that field of view given 
that it's in that field of view. 

There are four types of searches one might define as optimal. One 
might be interested in maximizing the total probability of detection when 
constrained to a given number of inspections (say K). If the incremental 
cost function Y(j»k) = c(j,k) - c(j,k-l) is a constant, then (after a 
suitable renormalization) one is demanding that P[f] be a maximum for C[f] 

< K. Such a search is termed totally optimal . If one demanded optimality 
for all K*l, 2, 3, . . . then the search is called uniformly optimal . A 
third type of search plan that one might consider is the search plan that 
maximizes the probability of detection with respect to the incremental cost 
and does so at every step of the search. Mathematically one finds the value 
of j which maximizes p(j)3(j,k)/Y(j,k) at each k. These searches are called 
locally optimal . Lastly one might entertain a search plan that minimized 
the total expected cost (i.e., was the fastest ) to find the target. 

The essential assumptions necessary to cast the asteroid search 
into the simplest form of the mathematical superstructure that Stone 
(1975) outlines are 

(1) That the asteroid is stationary (i.e., search rate high 
compared to the asteroid's angular speed), 

(2) That the search space is discrete (i.e., a fixed field of 
view) , 
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(3) That the allocation of effort is discrete (i.e., no 
favored fields of view), and 

(4) That y is aounded away from zero and p(j)b(J,k)/Y(.j,k) 

is a decreasing function of j (i.e., no free examinations 
of a field of view and the larger the search space the 
more difficult to detect). 

I do not believe that the physics or astronomy is strained by these 
strictures. In fact (5) y = constant is not unreasonable (i.e., the 
telescope moves smartly. The important point is that under these five 
limitations the totally optimal search plan, the uniformly optimal search 
plan, the locally optimal search plan, and the fastest searches are all 
Identical . Not only that, it can be explicitly exhibited. See Stone's 
text for the rigorous mathematical statements of the relevant theorems and 
their proofs. 


4. THE SEARCH PLAN 

I need just a bit more mathematics before I can exhibit the solution 
to the optimal search problem. The search plan C - ^ 2 * C 3 , ...) is a 

sequence of values e J for i = 1,2,3, ... These specify that the i'th 
examination be in field of view if the previous i -1 inspections failed to 
detect the asteroid in fields of view Let the set of all 

such search plans be denoted by H. Introduce the probability P[n, C] (and 
the cost Cfn, 5 }) of detecting the asteroid on or before the n'th examination 
while performing search plan f, e H (of the first n inspections). Finally, 
let r(j,n,C) be the number of scrutinizations out of the first n that are 
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placed in j'th field of view while following search plan C. A uniformly 
optimal search plan [for ~ Ij this is an unimportant normalization J 

I. H is one such that 

P[n,r*3 - max (P[n.r,] ; r h’l . n - 1, 2 K 

A locally optimal search plan is one such that is determined by 
[y ^ 0 necessarily] 

P(^i)3(€t. 1) _ P(j)e(j,l) 

Y(e,,i) y(j.i) 

and having determined the field of view for the first n-1 examinations 
(Ci* ^2**“^n-l^ field of view for the n'th one is determined from 

p(i)3(i,r(i>n-l,C) -H) _ p(j)6(j.r(j>n-l,g) +}) 

Y(i,r(i,n-l,C) + 1) y(J , r(j,n-l ,C) + 1 ) 

with Cp = i. Now define - r(Cj^, n, 5). The notation means that 
the n'th examination of the search plan i is placed in field of view 

and that it is the ^y, time that this field of view has been searched. 
The average cost to find the asteroid can be expressed in a variety of ways 
if the limit as n « of P[n,C] is unity; 

CO 

y(e) = E C[n,C] (P[n,d - P[n-l,d) 

n-1 

00 f, 

ci;!ai>CAL t:j 

OF POOR QUALITY 
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r.ince P[0,/J - (j. if y(j,k) " 1 then thir> rcducef; to 

P(0 - E d-PCn.f,]) 
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Now I can exhibit the solution explicitly. Under the assumptions 
outlined above if qj is the probability of detecting the asteroid after 
a single examination of field of view number j (given thau it is in 
field of view number j) then, as each inspection is an independent event, 
the incremental conditional ov'bability of detection 3(j.k) = b(j.k) - 
b(j,k-l) is given by 

B(j,k) = qj n-qj)'''’ for j e J , k = I, 2, . . . 

Normalize such that Y(j.k) = 1 V j e J, k = 1 , 2. . , . aqd suppose that an 
allocation f(j) has total cost (i.e., number of inspections) K. 

E f(j) = K 
jeJ 

The total probability of detection for this allocation of effort will be 

P[f] = E P(j)b(j,f(j)) = D p(j)[l-(l-q 
JeJ jeJ J 

Consider the search plan defined by: one makes the n'th inspection in 

field of view number i e J such that 
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/•\ /I n, C) fllOX p/ j 1 ; \ 

O-q,) "J,.J p(J)<, (, -q.) '■ 

J J 

Fhon f - f;* a„d K optimal (i„ .,|1 noptoc). Thin it d,„. 
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J' J qjd - qj) 


^ (j » If"! , (,) 
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5. SEARCH PLAN CONSTRUCTION 

I have already argued that the a priori target distriauU,, p(j) 
be approximated by a defective uniform distribution over the st irch e. 
In fact, p(J) V area (J)/4v. I have also argued tfv pe t , total 
cost function is homogenous over J and Indepeo, , ,.,a number of looks, 
Y(j,k) = 1 (in appropriate units). The probability of detection of the 
minor planet in field of view number j (given that the asteroid is there) 

Is qj. This depends principally on the apparent magnitude of the minor 
planet and the night sky background. Three effects tend to make asteroids 
fainter; atmospheric extinction, loss of brightness due to increasing phase 
angle, and increasing distance (heliocentric or geocentric). 

The extinction is modeled as usual, 


e = secz 


where t is the topocentric zenith distance and is the extinction per 
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unit air mass. I've used a value of 0.13 mag/air mass for e^. For the phase 
funct'ic'.i in magnitudes I've used the Gehrels and Tedesco (1979) results, 

B(1,0) = B(1,0) + 0.538 - 0.134 - 7^1, for \q\ < 7“ 

B(1,0) - B(1,0) - 0?. for |0l > 7" 

where B(1,0) is the absolute B magnitude and B(l,e) is the apparent magnitude 
corrected for phase angle 0. The parameter of the linear part of the phase 
function in magnitudes I = 0.039 mag/deg. The search plans in Figs. 2 and 3 
presume a geocentric distance of 0.5 A.U. at opposition. 

For asteroids very much brighter than our limiting magnitude (m - mj^ < -1 ) 
the probability of detection is essentially unity. For asteroids very much 
fainter than our limiting magnitude (m - m^ > l"*) the probability of detection 
is essentially zeroj see Fig. 1. Hence, the most interesting range from 
the point of view of planning a search is the regime jm - mj^| < 1/2'”. The 
search plans shown in Fig. 2 are for midnight on a winter solstice night 
and a B(1,0) = m^ - l"* (Fig. 2a), or m^ - 1/2'" (Fig. 2b). The search space J 
was chosen to be the 20^^ x 2*^ (declination x right ascension) area on the 
celestial sphere centered at opposition. Note that this prejudices the 
search plan towards the intuitively obvious region of the celestial sphere. 

(The latitude of our observatory is 33°49‘.) Each field of view of the search 
space is a square, two degrees per side and there are 150 fields of view in 
the search space. We look first in the field of view with the highest 
probability of detection and choose subsequent fields of view based on 
Eq. (1). A simple, repetitive enumeration through Eq. (1) determines the 
field of view order. Figure 3 shows the mj^ - 1/2'" case at midnight on a 
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summer solstice night. All of these illustrations were arbitrarily 
terminated as soon as each field of view of the search space had been 
examined (for clarity in presenting the plans diagrammatically) . 

A logical question to ask at this point is "How inefficient is 
the unplanned search relative to the optimal search?" For our searches 
we would've started at opposition and then spiralled outward until each 
field of view had beer examined once. When the asteroid is relatively 
bright (Figs. 2) this is roughly the same as the optimal search plan. 

The search plan exhibited in Fig. 2a has a cumulative probability of 
detection of 92.4%. For the search plan shown in Fig. 2b the cumulative 
probability of detection is lower, 85.3%, and the usual search plan is 
5.5% less efficient still. The comparison of the two search plans for the 
case of Fig, 3 is more complicated because we would've never repeated an 
examination of a field of view. With comparable effort to the optimal 
search plan, but randomly distributed over the search space, the optimal 
plan starts out more efficient and then becomes comparable to the repeated 
uniform in areal coverage one. This is typical of extended optimal search 
plans for medium bright objects. When one is trying to fully reach one's 
limits optimal search plans invest tremendous allocations of effort 
repeatedly near opposition (since each examination of a field of view is 
an independent event and the sought for asteroid is at the limits of 
detection). The plans tend to be factors of 2-4 times as efficient of 
the uniform ones. 
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6. GENERALIZATIONS 


It Is clear that any search for a fixed object, from variable stars 
to geosynchronous artificial satellites, can be cast into this formalism. 

It should be just as clear that this paper contains all of the essential 
mathematics for single searches of this type. Searches for moving objects 
and multiple observatory searches for the same moving objects can also be 
solved by similar methods. They are however much more difficult to formulate 
and specify especially since their theoretical structure is incomtlste. A 
more relevant problem is the multiple night search (by the same observatory) 
for a fixed object. One can plan such searches by an iterative algorithm 
that takes into account the (presumed) failure of the search. The posterior 
target distribution, given failure to detect, is updated by Bayes's formula 
and the conditional detection probability is appropriately modified too. 

In this fashion a whole week's worth of searching can be optimized. These 
techniques are applicable to all types of searches (x-ray bursters to comets), 
and can handle false targets, approximations to optimal plans by incremental 
means, etc. 
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Figure Captions 


Figure 1 


Figure 2 


Figure 3. 


. Probability of detection as a function of "distance" 
from limiting magnitude mj^. The functional form is 
constant/n + exp [5(m-mL)]}. 

. Search plans for (Fig. 2a) a bright and fainter (by oTS; 
Fig. 2b) asteroid at midnight on a winter solstice night. 
The number(s) in the boxes are which examinations of the 

optimal search plan this field of view of the search space 
was examined. 

Same format as Fig. 2 except at midnight on a sumner 
solstice night for m = m^^-1/2’". 


19 



-a/ 


ORIGHMAL I’Ain; Tj 
OF POOR QUALITY 


1.01 




NORTH 



RIGHT ASCENSION 



north 




ft 

\ 



o 

0 ) 


right ascension 




























¥ Sr/ 


ORlGiWft^ PAGE S3 

OF POOR QUALfTY 


H 

\ 


m 

D 


n 

m 

■ 


m 

mm 

■ 

« S 5 5 

t- (I •* 

f* 4- n r* 

4S 4n 

••' r» If 

}5j::}?2 

»- r» 4* 

4* 4*4 «» 

40 fs *- K 
1' » 40 

40 

^ o o p 
CO ( ) CO cv 

Cc c« «» 

40 p- ^ PC r. 

4 > P« f- r' <0 

r- t4 n n >t) 

i« 40 n o o 

<x CX 2 5 2 

o 4f> «» X» p- 

<0 XI in <0 x> 

(« rc PI «♦ XI 

4* P) X> P- 

3 « 4- « 
40 «0 XI XI 

<S « In o 
, *• <>♦ 

♦• « r> w 

<o a» f (.» 

♦- '») 40 

K «p 4f 5 

lO 4f »n 

4>» l*t 40 

lO 40 40 <»» 

40 »C 4 4 o 

«» <M t4) 

4T* .0 O 

4/1 </) w> »■ 
♦- rf> 

p C^ Sr» <n 

♦- (C 4f 

tO m 4» PC 
ttl 01 C' 40 

4* PC «» 

•» »» *0 P- 

PI PI IX *» 
P" PX <x «o 

PI PI m (»> p’ 
«• 40 C-. ,,| rx 
PC c« p» X» XI 

0 0 4» 

<'»)«)♦- 

<*» «0 •» 40 

2 5 n 5 

«• n tf) 

^ m 

<n 0 * m 

N O 40 
44| ^ O 
♦- n 

« CO «n 
<» ♦- «n 
*• m 

p- <n 40 

5522 

*- PC ^ 

2 5 2 5 

CX P* V 

rx 01 K m PI 

O ^ XT PC PC 

PC CX C*l ^ lO 

p* ri 10 01 X) 

^ p» P* 01 x» 

cn PI 40 V XI 

^ ♦ N 

•2S 

r» 8> 
rt #>♦ N 
»» rt 

^ 2 ^ 

04 40 o 

^ « 40 
*• ♦ 

»- »• lO 
PCO04 

#• nl ♦* 

^ Jf ® 

O CO PC 
40 9 PI 40 
PI 40 

01 40 PI 

19 O 40 40 
P* PC 40 

OAPCPip- 
40 PI cn f « 
P* PI PI 40 XI 

p* cn CX XI 40 

XI ^ ® K ^ 

CX PI PI ♦ XI 

rf> N K 

<n ^ 

^ n 

15 5 S 

4- ♦ 

2SS; 

*• n 

0) 40 00 
CO 40 
♦- 40 

40 40 (0 
A N 
*• 40 

<0 PI 4» 
p 40 K 
*• ♦ 

? ? 5 

p* PI 

01 cn ^ 

40 P» 01 01 

♦- CX ♦ 

m XI CO fs 4f 

o CX o ^ « 

^ IX PI 40 X) 

XI 40 CX 40 0 

® 0 • « PI 

CX cn rt ^ XI 

A « ^ 

^5!S 

40 40 
^ 40 40 

0* « o 

«> «0 40 
f 40 

^ «0 
fX 

<*> 

<S 

^ 40 


PI 40 XI 40 

40 40 6 o 

P* « X) 

8 rx o 5 5 
^ PI PI 9 XI 

♦» PC « 40 XI 

9 40 N XI CX 

PC CX cn ^ XI 

O N «0 


^ s s 

♦* W) 

<n 01 
cx 

PC 40 

3 5 

If N PX K 

^ 55 

CO CD p» 
PC CX 00 
p- PI 

88S2 

P4 PI XI 

33338 

CX PI ♦ XI 

PC PI X) CS XI 

« « rs XI N 

PC (X cn ^ XI 

m 


♦ <« 40 
*► 40 « 

«- 40 

O lO fH 
^ 00 40 
»» 40 

40 

525 

f- « 


5522 

40 40 01 40 m 

o CX g p* 40 

1 P* PC ^ 40 40 

g N C4| XI 4F 

# 0 ® • n 

PC « PI ♦ XI 

*28 


40 CO « 
M n o» 
^ n 

O 04 40 

? 

IX Q M 
PC 40 <X 

*• ♦ 

228 
p* PI 

CX ^ fs P* 
40 « PI 40 
♦P PI CO 

2 2 40 2 
P- PC 40 

2 « 40 o 

40 <n PI 4P XT 

p* PC PI ^ A 

8 2 2*2 
CX C«l ® ^ « 

S$8S 

n 40 

• 4> 40 


srs 

♦- n 

25 8 

^ n 

8582 

^ PI 40 

S588 

40 PI 4- P* 
01 ♦- 01 PI 
CX CX ♦ 

SS285 

IX CX PI ^ « 

(O PC PI A 40 
2 « 2 40 * 


c* « O fs 

•J2S2 

CO 01 n <n 

40 ♦ 4^ « 

40 IS fv M 

o K PC <n 

40 40 «» O 
p* 40 

2S|| 

*255 

♦P PC ♦ 

8833$ 

^ CX cn 40 5 

« 01 rs 01 CO 

3 40 p4 XI fx 

PC PC PI ♦ XI 

5S25 

PI # # A 

sag| 

CO 5J N ^ 

♦ « • r> 
4'* N N «0 

•• N ♦ 

•- PC 40 

SS58 

p* PC ♦ 

2 o S 2 

PC PC XT 

« CO 40 *- 01 

O CX ^ ^ 40 

*p PC PI ^ 40 

^ S t ♦: ® 

Px * PI ? S 

*: 45 « ^ ® 

3883S 

A A N IX 

3 A p* A 

4 « A A 

« 40 N 

Ot r« 40 N 
M N 4» 

«ssa 

N N «r 

4»C N ♦ 

S22S 

CN M 4» 

SS85R 

t- PC <n 40 « 

5s3l* 

NOO^t- 
O 40 40 PI N 
PX PC PI ♦ « 

2 * S 5 8 

PC PI PI ♦ XI 

2285 

PI PI t m 

U 

S§sH 

fs <n » o 40 
o rc f 40 

M 10 ^ 40 

<»» t- O C» M 

o m M o 40 

4- oc 40 ♦ 40 

^ Q ^ lO N 
O* 10 <0 40 

4^ t PC ca <0 

2 S 2 2 ‘S 

p* PC r» ♦ 40 

4 in N CO o 

PC in M 9 r> 

PC cx PI 40 40 

(Q 01 <n a CX 

xc CX CO 3 4«i 

PC PI PI ♦ tfl 

p> 40 m xf ^ 

SS2S5 

® 40 ♦- CJ 

2258 

¥ 

2 ^ 

8 1 

4a 

4 


0 } 

< 

UJ 


oc 

o 

z 


3 

O 


NoiivNuoaa 


Right ascension 




